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The results of the mean field studies on the effects of the backward scattering with large 
momentum transfer in a two-dimensional electron system are extended to the case with various 
types of the Fermi surface and coupling constants. It is found that the coexistent state of d- 
wave superconductivity, antiferromagnetism and 7r-triplet pair can be stabilized quite generally 
near half filling. 
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§1. Introduction 

In recent years the electronic states in two-dimensional systems in the copper oxide high-T c super- 
conductors have been studied intensively. Especially in the normal state near the optimal doping, 
a very flat dispersion of quasiparticle excitations around (n,0) and (0,7r), i.e., the extended saddle- 
point singurality, has been revealed in the angle resolved photoemission spectroscopy (ARPES) ex- 
periments.!' & This flatness of dispersion is more prominent than that obtained in the single-particle 
band calculations such as local-density approximation (LDA).i'& This behavior has been attributed 
to the many-body correlation effects based on the results of quantum Monte Cairo (QMC) sim- 
ulations,!^ and propagator-renormalized fluctuation-exchange (FLEX) approximation or to the 
dimple of the Cu02 planes from the recent LDA calculations .B It is the well-known fact that saddle 
points in the electronic structure result in the van Hove singuralities (vHs) in the electronic density 
of states (DOS). Especially for two-dimensional case, they produce logarithmic divergence in the 
DOS. Therefore, their existence near the Fermi level has important physical consequences. 

In the hole-doped case, such as YBa2Cu307_,5 (YBCO) or Bi2Sr2CaCu20s+5 (BSCCO), the Fermi 
level lies near (vr,0) and (0,7r).B This might be responsible for unconventional physical properties, 
such as the normal-state resistivity with T-linear temperature dependence^ or the superconduct- 
ing gap of the d x 2_ y 2 symmetry.il 0) On the other hand, in the underdoped region, a pseudogap 
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of the d x 2_ y 2 symmetry, which is consistent with the superconducting gap symmetry, has been 
observed around (ir,0) and (0,7r) above the superconducting critical temperature .0 It has 
been indicated that there is strong coupling between quasiparticle excitations near the flat band 
and collective excitations centered near (7r,7r).0 

In the electron-doped case, such as Nd2-xCe x Cu04 +< 5 (NCCO), however, the Fermi level lies far 
above (n,0) and (0,7r), and instead gets close to (±7r/2, ±7r/2). The Fermi surface agrees 

very well with LDA results & This might be related to conventional (or Fermi-liquid type) phys- 
ical properties, such as normal-state resistivity with T-square temperature dependence© or the 
superconducting gap which is nodeless on the Fermi surface (or BCS-like gap function)© 

Theoretically, various shapes of the Fermi surface can be reproduced by introducing not only 
nearest-neighbor hopping t but also next-nearest-neighbor hopping t'. Studies looking for the 
instabilities have been carried out with a special emphasis on (vr,0) and (0,7r) in the 2D Hubbard 
model, by use of the renormalization group methodHS and QMC© These studies have shown 
that d-wave superconductivity prevails over antiferromagnetism by the effect of t' . In our preceding 
letter,© we have studied possible ordered states with a special emphasis on the effects of the 
backward scattering processes with large momentum transfer between two electrons near (tt,0) and 
(0,7r) by introducing g\ and g% processes, which correspond to 'exchange' and 'Umklapp' processes, 
respectively. The mean field phase diagram for a special case, g\ =53, and for the YBCO type 
Fermi surface, was determined. In this paper, results of more detailed studies for general cases with 
various shapes of the Fermi surface and for several choices of g\ and 53 will be presented. In §2, 
we introduce the model Hamiltonian and three types of order parameters, i.e., ci-wave Cooper pair, 
Neel order and 7r-triplet pair. In §3, we determine the phase diagram in the plane of temperature, 
T, and the hole (or electron) doping rate, 5 (or x). §4 is devoted to conclusion and discussion. 

We take unit of h = k& = 1. 

§2. Model 

We consider a two-dimensional square lattice with the kinetic energy given by, 



<1]>(T 



per 

where tij is the transfer integral, Ci CT (c| ) is the annihilation (creation) operator for the electron on 
the i-th site with spin a, and [i is the chemical potential. The energy dispersion £ p = e p — [i is given 
by 

e p = —2t(cosp x + cos p y ) — At 1 cosp x cosp y 

-2t"(cos 2p x + cos 2p y ), (2.2) 
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including t (nearest neighbor), t' (next nearest neighbor) and t" (third neighbor), as shown in 
Fig. 1. We take t as the energy unit, i.e., t = 1 and the lattice constant is also taken as unity. The 
energy band e p has saddle points at (±7r,0) and (0,±7r). The number of independent saddle points 
in the 1st Brillouin zone is equal to 2. 
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Fig. 1. The transfer integrals on a 2D square lattice; t (nearest neighbor), t' (next nearest neighbor) and t" (third 
neighbor) . 



The 1st Brillouin zone is usually taken as shown in Fig. 2 (a). In this case, however, the saddle 
points lie on the zone boundaries. In order to treat the scattering processes between two electrons 
around the saddle points unambiguously, it is convenient to choose the 1st Brillouin zone inside 
which the saddle points lie, as shown in Fig. 2 (b). It consists of the region A and B including 
Qa = (tt,0) and Qb = (0,vr), respectively. We note that w p = sgn(cosp y — cosp x ) = +1(— 1) for 
p GA(B), which determines the symmetry of the superconducting order parameter, as we shall see 
later. 

Now we introduce effective interaction between two electrons in the region A and B, just as in 
the g-ology in one-dimensional electron system. Here we treat only the backward scattering with 
large momentum transfer, shown in Fig. 3, i.e., g\ and 53, in analogy with normal and Umklapp 
processes in our previous work for one-dimensional electron system.!!' 1 We take only interactions 
between two electrons with antiparallel spins into account, as in the Hubbard model, i.e., g\ = g%± 
and 53 = g3±. We treat the above Hamiltonian in the mean field approximation. We are interested 
in the repulsive case, 

< gi < 53- (2-3) 



We consider three types of the scattering channels: 
(1) Cooper-pair channel 
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Fig. 2. The 1st Brillouin zone, (a) original one. (b) our choice. 
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Fig. 3. The backward scattering processes with large momentum transfer. The solid and dashed lines stand for 
electrons in the region A and B, respectively. 



This channel consists of the scattering processes as shown in Fig. 4 (a), where a pair of two 
electrons with total momentum 2Qa is scattered to that with total momentum 2Qb and vice 
versa. This is included only in g% processes. We note that this channel is represented in the 
original 1st Brillouin zone as shown in Fig. 4 (b), where there exist only normal processes between 
two electrons with total momentum equal to zero. Therefore, we naturally introduce the order 
parameters of the Cooper-pair, 

A A = E'fc < c QA-HCQ A+ k-\ >, 

(2.4) 



\k x \ + \k y \<7r 



k , , -|fe|+|*y|<ir (2vr) 5 

It is to be noted that Qa ± k and Qb i k always lie in the region A and B, respectively, for k 
satisfying \k x \ + \ k y \ < tt shown in Fig. 5. This is the reason for the choice of the Brillouin zone as 
shown in Fig. 2 (b). 
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Fig. 5. The region satisfying \k x \ + \k y \ < tt (shaded area). 

For <?3 > 0, the order parameters, A^ and Ab, can be non-zero only in the case, 

Ai = g 3 A A = -g 3 A Bl (2.5) 
i.e., the superconducting gap has d x 2_ y 2 symmetry: 

A X 2_„2 = W P < c -vi c v\ >> ( 2 - 6 ) 
v 

= A A - Ab, 
= W 

93 
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We note that this superconducting gap function is nodeless in spite of d x z_ y i symmetry, i.e., 
only its sign changes across the boundary between the region A and B and its magnitude is 
constant all over the 1st Brillouin zone, which means that it changes dis continuously along the 
region boundary. This can be more clearly seen from the fact that we can rewrite ( |2.5| ) as 
A(p) = g^^-iWp = <73Aisgn(cospy — cosp x ). It is due to our choice of the 1st Brillouin zone 
as Fig. 2 (b) and introduction of the 'Umklapp' scattering g% with momentum dependence ignored. 

This channel is the same as 'pair-tunneling' of two electrons located around Qa and Qb, which 
favors d x 2_ y 2 pairing, in the 2D Hubbard model.!!' These interaction processes are not present 
as important factors in the recent model by Assaad et ah, who have introduced the additional 
interaction expressed as the square of the single-particle nearest-neighbor hopping.© 
(2) density-wave channel 

This channel, consisting of the scattering processes with momentum transfer equal to Q = (ir, ir), 
as shown in Fig. 6 (a), is included in both g\ and 33 processes. In the original 1st Brillouin zone, 
this is represented as shown in Fig. 6 (b), where there exist both normal and Umklapp processes. 
We introduce staggered carrier density of each spin, 
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Fig. 6. The density-wave channel (a) in our Brillouin zone and (b) in the original Brillouin zone. 



A T = T!k < C Q B +^ C QA+k*\ >, 

< 

t A i = E' k < c q b -u c Qa-h > ■ 

For (|2.3|), they can be non-zero only in the case, 



(2.7) 



Af = —A* = real, 



A 2 = 5sA T = -3 s A i; 



(2i 



where g s = g\ + g%. This is the commensurate spin-density- wave (C-SDW) or antiferromagnetic 
(AF) state. The staggered magnetization (Neel order) M is given by 

M = ^]T(-ir<Sf >, (2.9) 

i 

\ ^ t t 

~ 2 ^ < c p+QT c pT ~~ c p +QI c pI >' 
p 

= Re(A, - A,), 

= 1a, 

In this paper, we consider the commensurate case only. The validity of this assumption will be 
discussed in §4. 
(3) 7r-pair channel 

This channel, consisting of the scattering processes between two electrons with total momentum 
equal to Q, as shown in Fig. 7 (a), is included only in g\ processes. In the original 1st Brillouin zone, 
this is represented as shown in Fig. 7 (b), where there exist both normal and Umklapp processes 
(only the latter case is shown in the figure). We introduce the order parameters of 7r-pair, 
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Fig. 7. The 7r-pair channel (a) in our Brillouin zone and (b) in the original Brillouin zone. 



A + = E'fe < CQ A -klCQ B+k ^ >, 

A - = E'fc < c Q s -fc|CQ A+fcT > . 
For g\ > 0, they can be non-zero only in the case 

A 3 = #iA+ = -#iA_, 



(2.10) 



(2.11) 
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i.e., n-triplet pair exists, 



^triplet 



< c_ 



(2.12) 



pa 



= -2(A + -A_), 
= -±A 3 . 

It is necessary to decouple the Hamiltonian by use of this order parameter because if dSC and AF 
coexist 7r-triplet pair always results. This important fact has been indicated by Psaltakis et al. in 
slightly different context. This 7r-triplet pair was not taken into account in the calculation based 
on the slave-boson mean field approximation in the 2D t-J model,@) where only dSC and AF are 
simultaneously taken into account. 
By use of the above three order parameters, we obtain the mean field Hamiltonian as follows: 

h mf = Y!4 M kA+E c , 



~ ( C QA+kV c QA-k[i c Q s +fc-p c Qs-fc|), 



M, 
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There are four quasiparticle energy bands, ±E + (k) and ±E_(k), 

E±{k) = J + l A il 2 + ( A 2) 2 + I A 3 | 2 ± A(k), 



(2.13a) 
(2.13b) 

(2.13c) 

(2.13d) 
(2.13e) 

(2.14a) 



A(k) 



\ 



(afc - b k f 



a k + b k 



+ |A 3 



+ [(o fc + 6jfc)(A2)-2Re(A|A 3 )] 1 



(2.14b) 



where we note that ak <-> bk-, A(k) <-> A(k) and E±(k) <-> E±(k) as we exchange k x and k y . 
The self-consistent equations are given by 

Ai = ^V'y ^-tanh^l 
4 ^ £±E a {k) 2T 



A, 



A, 



tE'e-^^^Ia 



|Ai + a^- [Re(A*A 3 ) - a fc A 2 ] j , 
(a k -b k ) 2 



k a=± 



E a {k) 



2T 



tanh ■ 



2T 



A, 



n 



2T 



1 — a 
1 + a 
l + a 



2A(k) 

{ak ~ b k ) 2 
2A(k) 

{ak ~ b k ) 2 
2A(k) 



(2.15a) 



- «|^y [Re(A*A 3 ) - a fe A 2 ]| , (2.15b) 

2A 1 
+ a-^y [Re(A^ A 3 ) - a fe A 2 ] I , (2.15c) 

OA ~) 

- [Re(A*A 3 ) - a fc A 2 ] \ ,(2.15d) 
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where n is the electron filling, which is related to the hole doping 5 and the electron doping x by 
5 = 1 — n and x = n — 1, respectively. These equations satisfy J^- = = = 0, where F is 
the free energy per one lattice site in the canonical ensemble given by 



F = -2T^' {log [ 2 cosh ^ 
+E C + fin. 

The electronic density of states p{e) is given by 



+ log 



2 cosh 



E-(k) 



2T 



(2.16) 



P(<0=E' E S(e-a'E a (k)) 

£«(*:) t°* + ° 2A(fe) 



1 + 



OA 

_ [Re(A^A 3 ) - a fc A 2j 



(2.17) 



which is related with n by n = J^° 00 dep(e)f(e,T), where f(e,T) = (e e/,T + 1) 1 is the Fermi 
distribution function. 

§3. Results 

The self-consistent equations (2.15) are solved numerically and the phase diagram has been 
determined in the plane of temperature T and the doping rate 5 or x. It is easily understood from 
(2.15) that [1] if two of the three order parameters coexist, another one always results and [2] for 
gi/t = 0, 7r-triplet pair cannot exist (A3 = 0), i.e., there exists no coexistent state and dSC prevails 
over AF. We note that it is energetically favorable to take both Ai and A3 to be real for the 
coexistent state. We choose t' and t" so that the bare band width W is equal to 8t and reproduce 
various shapes of the Fermi surface. We fix the value g^/t = 5.0. 

In the numerical calculation of the DOS (2.17), we approximate 5-function by 

*(e)=- hm |-/(e,77)= hm {4r/cosh 2 f }-\ (3.1) 

with a finite value n/t = 0.01. 

First we include t only (t'/t = t"/t = 0). In this case, there is perfect nesting property e p = — e p +Q 
or bk = —cik — 2/Lt. Since the quasiparticle dispersion depends on k only through a^, it is convenient 
in the numerical calculation to rewrite (2.15), (2.16) and (2.17) by use of the following identity for 
any function g of at , 

E' ' 9(ak) = g / de Po(e)5(e ~ f 1 )- ( 3 - 2 ) 

k 

Po(e) is the bare DOS given by 



Po(e) = 2<5(e-C P ), 
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where K is the complete elliptic integral of the 1st kind K (k) = Jq —?= ^ — ^— and 9 is the 

y(l— x )(1— k x ) 

step function. 

We show the bare Fermi surface in Fig. 8 for a few values of the hole doping 5 (it is sufficient 
to consider the hole-doped case due to the particle- hole symmetry). The saddle points lie on the 
Fermi surface in the half- filled case. 



t'/t=t"/t=0 




D TC 



Fig. 8. The Fermi surface in the case without t' and t" for a few values of the hole doping 8. 

We show the phase diagram in Fig. 9 for two different choices of the coupling constants, g\ = 53 
and g\ = In the half-filled case, as the temperature is lowered, AF gets stabilized, and since 

there is no free carriers in the presence of the SDW gap, dSC cannot arise even in T = 0. On 
the other hand, near half filling, there are free carriers even in the AF state, and dSC also can 
be stabilized in the low temperature region. This result is similar to the spin-bag picture© It 
is important that the coexistence of AF and dSC always results in -zr-triplet pair. Such a close 
relationship between dSC and AF has been indicated by SO(5) theory© For smaller g\, AF region 
and therefore the coexistent region are reduced. 

We note that for gi = and 5 = 0, due to the perfect nesting property, Ai and A2 are determined 
by the same equation, i.e., AF and dSC are degenerate. 

In Fig. 10 we show the DOS in the coexistent state near the optimal doping which gives the 
highest onset temperature of dSC. There exists energy gap in the DOS at the Fermi level. 

Next we consider the hole-doped case with the YBCO type Fermi surface. We choose t'/t = — 1/5 
and t" /t = 1/7. The Fermi surface is shown in Fig. 11 for a few values of the hole doping 5. In this 
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Fig. 9. The phase diagram in the case without t' and t" for gzjt = 5.0 and g\/t = 5.0, 2.5. 
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Fig. 10. The electronic DOS in the case without t' and t" in the coexistent state for gs/t = gi/t = 5.0 and T = 
near the optimal doping 5 = 0.21. The Fermi level lies at e = 0. 



case, the saddle points lie on the Fermi surface near 5 = 0.44. 

We show the phase diagram in Fig. 12 for two different choices of the coupling constant, g\ = 
and <7i = gs/2. For 5 = 0, only AF can be stabilized, and near half filling there exists coexistent 
state of dSC, AF and -/r-triplet pair. For g\ = g^, compared with the perfect nesting case, the effect 
of nesting near 5 = is suppressed and the onset temperature of AF is lower. On the other hand, 
near 5 ~ 0.4, the saddle points lying near the Fermi surface lead to the enhancement of the order 
parameters of both dSC and AF. 
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Fig. 11. The Fermi surface for the YBCO case for a few values of the hole doping 8. 




Fig. 12. The phase diagram for the YBCO case for gz/t = 5.0 and gi/t — 5.0, 2.5. 



In Fig. 13 we show the DOS in the coexistent state near the optimal doping. The DOS has an 
energy gap at the Fermi level. 

Finally we consider the electron-doped case with the NCCO type Fermi surface. We choose 
t'/t = —1/3 and t" /t = 1/10. The Fermi surface is shown in Fig. 14 for a few values of the electron 
doping x. In this case, the saddle points lie on the Fermi surface near x = 0.28. 

We show the phase diagram in Fig. 15 for two different choices of the coupling constant, g\ = g% 
and g\ = g%/2. In this case the saddle points are not located near the Fermi surface and there is a 
suppression of the onset temperature of AF and dSC compared with the above two cases. 
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Fig. 13. The electronic DOS for the YBCO case in the coexistent state for gz/t = gi/t = 5.0 and T — near the 
optimal doping 8 = 0.30. The Fermi level lies at e = 0. 




Fig. 14. The Fermi surface for the NCCO case for a few values of the electron doping x. 



In Fig. 16 we show the DOS in the coexisting state near the optimal doping. The energy gap lies 
at the Fermi level in the DOS. 

§4. Conclusion 

We have studied possible ordered states of interacting electrons on a square lattice with a spe- 
cial emphasis on the backward scattering, g\ and g%, i.e., 'exchange' and 'Umklapp' processes, 
respectively, between two electrons around (vr,0) and (0,7r) for various shapes of the Fermi surface. 

We focus on the case that gs, which is related with the superconducting state with d x i_ y i type 
of the gap, is about half of the bare band width, W/t = 8.0, i.e., (/3/t=5.0. For g\ = 0, d-wave 
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Fig. 15. The phase diagram for the NCCO case for gz/t = 5.0 and g\/t = 5.0, 2.5. 
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Fig. 16. The electronic DOS for the NCCO case in the coexistent state for gz/t = gi/t = 5.0 and T — near the 
optimal doping x = 0.22. The Fermi level lies at e — 0. 



superconducting state (dSC) competes with and prevails over antiferromagnetism (AF) even in the 
half-filed case. For < g\ < g% with gi/gs not so small, only AF can be stabilized in the half-filled 
case and coexist with dSC near half filling. Especially in the coexistent state, 7r-triplet pair always 
results and the DOS has an energy gap at the Fermi level, independent of the hole or electron 
doping rate. The above conclusion is qualitatively independent of the shapes of the Fermi surface. 

As we have indicated in §2, we have assumed C-SDW throughout this paper. In general, however, 
away from half filling, the incommensurate SDW (IC-SDW) should be expected with Q / (ir, tt). 
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Therefore, it is not clear whether this coexistent state can be stabilized or not when we include IC- 
SDW. In fact, the instability towards incommensurate IC-SDW depends strongly on momentum- 
dependence of both the zero-frequency susceptibility, xo(o)i an d effective interaction. We note, 
however, the following facts: [1] There exists the region where Xo(l) shows maximum at q = (vr,7r) 
near half filling at low temperature, in the present choices of the Fermi surface. [2] g\ and g% 
can be considered to have maximum at q = (it, it) (q stands for the momentum transfer of one 
electron), although the g-dependence of g\ and 53, which is not explicitly indicated in this paper, 
is irrelevant under the assumption of C-SDW of interest. These facts result in RPA susceptibility 
showing maximum at q = (n, 7r) leading to C-SDW at least near half filling at low temperature 
where the coexistent state is possible. 

With respect to the superconducting gap symmetry, in the case that (±7r/2,±7r/2) lies near the 
Fermi surface, not only d x 2_ y 2 component but also other ones can mix. Recent QMC calculations 
show that not only d x 2_ y 2 but also d xy pairing correlation are enhanced in the 2D Hubbard model,© 
and new ordered state with a mixing of extended s-wave component with total momentum (vr,0) and 
(0,7r) and the usual d x 2_ y 2 component has been discovered in the 2D t-J model, in the mean field 
and Gutzwiller approximations.© Their results are closely related to the fact that the original 
d x 2_ y 2-w&ve superconducting gap has nodes around (±7r/2,±7r/2), i.e., the fact that Umklapp 
scattering between two quasiparticles on the Fermi surface around these points is possible may lead 
to form the full gap all over the Fermi surface. Our d x 2_ y 2 -wave superconducting gap, however, has 
no node, i.e., changes discontinuously along the boundary between the region A and B, as already 
discussed in §2. Therefore, in our model, we cannot discuss the mixing of another component in 
the superconducting gap function in addition to the d x 2_ y 2-wave component. This is the future 
problem. 
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